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Abstract 

We briefly introduce the generic framework of Disordered Elastic Systems (DES), giving a short 'recipe' of a DES 
modeling and presenting the quantities of interest in order to probe the static and dynamical disorder-induced properties 
(S| of such systems. We then focus on a particular low-dimensional DES, namely the one-dimensional interface in short- 
^ ranged elasticity and short-ranged quenched disorder. Illustrating different elements given in the introductory sections, 
Q we discuss specifically the consequences of the interplay between a finite temperature T > and a finite interface width 
^ > on the static geometrical fluctuations at different lengthscales, and the implications on the quasistatic dynamics. 

^-H Keywords: disordered elastic systems, interfaces, glassy phenomena 
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Could some features of experimental systems as dissim- 
ilar at a microscopic level as superconductors, magnets, 
ferroelectrics, fluids, paper, or two-dimensional electron 
gases, be described by the same equations at a macroscopic 
level? All those systems may actually display emergent 
structures such as interfaces (e.g. ferroelectric [1, 2, 3] or 
ferromagnetic [4, 5, 6] domain walls, contact line in wet- 
ting experiments [7] or propagating cracks in paper and 
thin materials [8]) or periodic systems (typically vortex lat- 
tices in type-E superconductors [9], classical [10] or quan- 
tum [11] Wigner crystals, or electronic crystals displaying 
charge or spin density waves [12, 13]). 

One can either describe them using ah initio predictions 
combined to a Landau approach, where two phases com- 
pete with each other at their common boundary (the com- 
plexity of a numerical approach increasing considerably 
with the system size), or rather take a radically opposite 
point of view by skipping the specific microphysics and 
focusing exclusively on the boundary, defined by the shift 
of the order parameter. Such an emergent structure can 
then be described as a fluctuating manifold or periodic sys- 
tem supported by a disordered underlying medium, in the 
generic framework of disordered elastic systems (DES). 

Thereafter we recall briefly the basic features of DES by 
giving first a short recipe of a DES model based on two 
competing physical ingredients: elasticity and disorder, 
blurred by thermal and/or quantum fluctuations. Then 
we list the main observables of interest in order to probe 
the disorder-induced metastability present in those sys- 
tems, and to address the two main questions which arise 
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regarding their resulting glassy properties: what can we 
learn by the study of its statics versus its dynamics, first 
via the characterization of its geometrical fluctuations and 
secondly via its response to an external force? 

Indeed, from the 1970s' and Larkin's work [14], we know 
that there could not exist a perfectly ordered solid in pres- 
ence of disorder, so how does the addition of disorder 
change the nature of a pure system? Finally we focus on 
a particular low-dimensional DES and study the static ge- 
ometrical fluctuations of a one-dimensional interface, via 
an analysis of the interplay between thermal fluctuations 
and a finite width of the interface in its roughness. 

Those short notes are not meant to be exhaustive, but 
rather to give a pedagogical and somehow practical intro- 
duction to the field, aimed at theoreticians but also at ex- 
perimentalists who might be interested in DES modeling. 
We focus essentially on the case of interfaces^ but most 
concepts can be extended to periodic systems^ and more 
details and references can be found for example starting 
from the existing reviews [13, 15, 16]. 

2. DES modelling: a recipe 

In the generic framework of DES, very few physical in- 
gredients are required in a minimal version of such a model. 
Thereafter we briefly sketch their concrete implementation 
for interfaces, but those considerations remain valid for pe- 
riodic systems [17]. 

2.1. Dimensionality and class of DES 

First of all one has to identify the dimensionality of 
the system {d being the internal dimension of the system, 
m the number of its transverse components, and D the 
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dimension of its embedding physical space) and whether 
it is a manifold {d-\-m = D with d = D — 1 for interfaces) 
or a periodic system. 

For example, a ID interface and a single vortex 
are manifolds respectively with {d = m = 1^ D = 2) and 
{d = l,m = 2,D = 3), whereas an Abrikosov vortex lat- 
tice and a 3D Wigner crystal are periodic systems with 
{d = D = 3,m = 2) said {d = m = D = 3). Note that the 
dimensionality and the class of DES might actually change 
with respect to some parameters of the system, as it can 
been addressed experimentally e.g. for ferromagnetic do- 
main walls (ID to 2D interfaces crossover) [18] or for vor- 
tices in superconductors (vortex lattice to individual vor- 
tices). 

2.2. Physical space of coordinates (z^x) 

The description of the physical space embedding the 
DES can then be split in two sets of coordinates: z G 
denotes the internal coordinates of the system (e.g. the 
position along a polymer or a given point in a lattice), and 
X eVx its transverse coordinates. 

For an analytical treatment of interfaces we typically 
assume that they live in an infinite and continuous phys- 
ical space so x Vx is taken as x R^. However, a 
physical or numerical realization of such a system is al- 
ways supported by a microscopically discrete sublattice of 
parameter 1/A, ultimately the crystal in a solid, and more- 
over lives in a finite box of typical size L with boundary 
conditions which have to be defined (possibly periodic or 
free). So in the comparison between analytical predictions 
and experimental or numerical results, corrections due to 
finite size effects and to the translation from the discrete 
to the continuous limits are a priori expected. 

Once the disorder is averaged out, a translational space- 
invariance is recovered, suggesting a description in Fourier 
space q with an ultra-violet 1/A and an infra-red 1/L cut- 
offs, which are always present in physical DES. However, 
from an analytical point of view, they are either irrelevant 
and thus skipped, or they are conveniently reintroduced in 
order to cure non-physical divergences in computations. 

2.3. Univalued displacement field 

In the absence of disorder, an elastic system tends to 
minimize its distortions, thus would typically be fiat for 
an interface or characterized by a single reciprocal vector 
for periodic systems. A given configuration of a DES is 
characterized by a univalued displacement field G 
with respect to such an equilibrium configuration (e.g. fiat 
or periodic) of the pure system as illustrated in Fig.l. 

The definition of this reference configuration is actually 
crucial but potentially tricky in experiments, in particular 
in certain cases where the equilibrium configuration is not 
a straight line. For example, how to define unambiguously 
the center and the mean radius of a 'dotted' ferromagnetic 
domain if it is far from being perfectly circular? Defects 
such as overhangs and bubbles for interfaces, or topological 




Figure 1: Definition of the displacement field for a ID interface 
{d = m = 1) superimposed over its surrounding smooth random po- 
tential V{z^ x) (in a weak disorder limit), is univalued only in the 
absence of bubbles or overhangs for interfaces. 

defects in periodic systems are still missing in the above 
DES description, since they hinder the definition of a uni- 
valued displacement field u^^ at the core of the concrete 
implementation of statistical averages both thermal and 
over disorder. 

2.4- Elasticity 

The elastic energetic cost of distortions iV ^u^) is given 
by the elastic Hamiltonian Hei [u]. For an interface, if the 
elasticity is short-ranged^ it is essentially proportional to 
the length of the interface Hei [u] (x J d^z ^/Y^^~{y^^^ 
and in the elastic limit of small distortions HV^il^ll <C 1 
it reduces to the quadratic form (thus analytically user- 
friendly) 

HAu] = lj^d''z-{V^sf (1) 

with the elastic constant c being the elastic energy per unit 
of length, and that corresponds in reciprocal space to an 
elastic energy cq^ per Fourier mode u^. 

If the elastic limit is broken, additional terms are a priori 
needed in the perturbation expansion V^TTTVJ^. Nev- 
ertheless the elasticity can also be effectively long-ranged^ 
as it is the case e.g. for contact lines in wetting experi- 
ments because of the fiuid surface tension [19, 20]. The 
effective energy per Fourier mode is then c\q\^ with /i 7^ 2, 
or the elastic Hamiltonian can include phenomenological 
additional terms as e.g. for the description of the out-of- 
equilibrium depinning regime (cf. section 3.2) [21, 22]. 

2.5. Disorder 

The 'disorder' accounts for the effects of the inhomo- 
geneities inherent to any real physical medium, and it is 
defined as a stochastic variable V{z., x) with a given sta- 
tistical distribution V [V] and the ad hoc disorder average 



2 



(denoted O for an observable O). It adds a random part to 
the DES Hamiltonian Hbes = '^ei+'^dis with the disorder 
Hamiltonian 

n^islu^V] = [ d^z-V{z,u,) (2) 

First, if this stochastic variable has a dynamics much 
slower than the dynamics of the DES, the disorder is 
quenched (as e.g. atomic terraces on thin epitaxial films) 
otherwise it is annealed (as e.g. itinerant oxygen vacancies 
in superconductors). That distinction essentially imposes 
the sequence of disorder and thermal averaging of observ- 
ables in computations. 

Secondly, disorder can either be dominated by a few in- 
dividual pinning centers in the strong disorder limit, or 
by many weak impurities in the weak disorder limit. In 
the latter case, the collective behavior of the impurities 
conspire (by the central limit theorem) to give rise to a 
smooth random potential F(£, x), as illustrated in Fig.l, 
with a Gaussian distribution V [V] . Note however that 
this 'weak disorder' limit can be realized in experimen- 
tal DES only because the disorder can actually vary on 
characteristic lengthscales which are much smaller than 
the lengthscales of interest of the DES description such as 
the lattice spacing of periodic systems (as e.g. for the 3D 
macroscopic Wigner crystal in colloids [23]). 

V [V] is thus fully characterized by its two first cumu- 
lants, namely its mean value V{z^ x) = and its variance 

V{z,x)V{z',x') = R^^{z-z')- R^^ [x-x') (3) 

The disorder is usually assumed to be uncorrelated along 
the internal direction z oi the interface with R^^ {z) = 5{z)^ 
whereas the disorder correlator R^^ {x) is the quantity ex- 
tensively studied in the Functional Renormalization Group 
(FRG) approach of DES (cf. section 3.3). If this correla- 
tor decreases sufficiently fast, then the disorder is short- 
ranged or random-bond (RB: the interface couples only 
locally to the surrounding disorder), else it is long-ranged 
or random-field (RF: the system is sensitive to the disor- 
der in all the physical space). In the section 4 we examine 
specifically what are the consequences of a finite RB dis- 
order correlation length > on the static properties of 
a ID interface. Of course in an experimental system all 
types of disorder could be present, but the dominant dis- 
order defines its universality class. 

2. 6. Internal structure of DES 

At last, the internal structure of physical DES can be 
crucially relevant for its properties, and it is actually al- 
ready partially encoded in the disorder correlator and its 
correlation lengths and (^x- 

The implications of an internal degree of freedom^ as 
e.g. the phase shift in Neel versus Bloch domain walls, are 
challenging theoretically, but such an additional physical 
ingredient in DES modelling might be needed to account 
for the physics of experimental systems, as it has been 



recently addressed in the context of spintronics' nanowires 
[24]. 

3. Observables as probe of disorder: statics versus 
dynamics 

There are essentially two basic questions that can be 
addressed regarding DES and their experimental realiza- 
tions: how do they look like and how do they respond when 
one pulls on them. By comparing the theoretical and nu- 
merical predictions of generic DES to measurements on 
experimental setups, it is possible on one hand to test the 
adequacy of a DES modelling, and on the other hand to 
identify the universality class of dimensionality, elasticity 
and disorder of a particular physical realization, and then 
to extrapolate to its other possible disorder-conditioned 
features. 

From an analytical point of view, the competition be- 
tween elasticity, disorder and thermal/quantum fluctua- 
tions is treated via two statistical averages of an observable 
O, respectively the thermal/quantum average (O) and the 
disorder average O. The two underlying assumptions to 
this procedure are the ergodicity and a self- averaging dis- 
order^ which imply essentially that the computed quantity 
(O) should match the measurement of the observable O on 
an equilibrated and sufficiently large experimental sample 
(here with quenched disorder). 

Thereafter we discuss briefly the statics and the dynam- 
ics of DES, and then we mention the main methods used to 
tackle them. For short pedagogical reviews on the subject 
and further references, see [15, 16]. 

3.1. Statics: geometrical fluctuations and roughness 

In statics the main information that can be accessed ex- 
perimentally is the configuration of the system, described 
by the displacement field defined with respect to a given 
reference configuration (which is somehow arbitrarily cho- 
sen) as in Fig.l. 

The geometrical fluctuations of can be quantitatively 
characterized by the probability distribution function 
(PDF) of the relative displacements Au^{r)= u^-^f — at 
a given lengthscale r. Assuming rightfully a translational 
invariance after disorder-averaging, this PDF V{Au{r)) 
can fairly enough be approximated by a widening Gaussian 
[25] as illustrated in Fig. 2, and its main feature is thus its 
variance, namely the roughness function B{f) = (A'u(r)^) 
or its corresponding structure factor S{q): 

Bi^= f j^2[l-cos{q-r)]S{q) (4) 

Those quantities are thus the mere two-points correlation 
functions of and its Fourier transform, for which we are 
hopefully well-equipped for an analytical treatment. More- 
over, if the system displays a scale invariance on a whole 
range of lengthscales, it is expected to behave logarithmi- 
cally or rather as a power law B{r) ~ r'^^ or S{q) ~ q~^~'^^ 
defining the roughness exponent (. 
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Figure 2: Definition of the relative displacements Aii^(r) for a ID 
interface and their translational-invariant PDF V{Au(r)) after the 
disorder average. If V{Au(r)) is Gaussian, its main feature is its vari- 
ance {Aii(r)2) = B(f), namely the interface roughness as a function 
of the lengthscale. 



For periodic systems, the displacement field is again de- 
fined with respect to a reference 'equilibrium' configura- 
tion, this time a lattice of reciprocal vector Kq. The anal- 
ogous quantity to the roughness is then the translational 
order correlation function (r) defined as: 



(5) 



and its Fourier transform the structure factor 5'(^o + q). 
If the theory is close to a Gaussian as in [17], it is actu- 
ally simply related to the roughness since we have then 

C^^(f) = exp[-f B(r)". 

In practice, in order to measure the roughness and its 
related quantities, one has first to take a picture (a 'snap- 
shot') of the system, then chose a reference configura- 
tion (typically flat for an interface or a perfect lattice for 
periodic systems) and define the corresponding displace- 
ment field u^^ and finally compute the appropriate func- 
tions 5(r), C{r) or S{q). In a logarithmic representa- 
tion, the roughness usually exhibits a power-law (or loga- 
rithmic) behavior B{r) ~ A{c^ T, ^) • r^*^ up to a certain 
lengthscale above which it saturates; the main focus then 
is on the value(s) of the corresponding roughness expo- 
nent (s) as a signature of dominant physics at a given 
lengthscale range depending on the universality class to 
which the system belongs. There are however other im- 
portant features also at our disposal to probe the disorder- 
conditioned properties of the system, such as the power- 
laws prefactors and their possible temperature-dependence 
e.g. A{c,D,T,() - T^^ (with the thorn exponent b), the 
crossover lengthscales including the saturation lengthscale 
itself and last but not least the possible non-Gaussianity 
of the PDF V{Au{f)) (that could then be due to finite- 
statistics artefacts or to a physical origin) [8] . The special 



case of a ID interface with a short-ranged elasticity and a 
RB disorder is discussed at length in the section 4. 

Such an analysis is of course also possible in order to 
study the dynamical geometrical fluctuations of the sys- 
tem, but that requires a time-resolved measurement tech- 
nique, and the possibility to take an actual snapshot of 
the system, which is not always the case (e.g. for surface- 
scanning techniques such as the Scanning Tunnelling Mi- 
croscopy or Atomic Force Microscopy). 

3.2. Dynamics: velocity-force characteristic 

The knowledge of the statics of a disordered system is 
not a sufficient criterion to decide on its glassiness, since 
it may look quasi-ordered despite the disorder (as e.g. pe- 
riodic DES displaying power-law decaying Bragg peaks). 
Because of its disorder-induced metastability, it may more- 
over be difficult to even discriminate equilibrium from out- 
of-equilibrium configurations, so it is necessary to study 
specifically the dynamics of the system. 

The characterization of the response of a disordered sys- 
tem under an external force is fundamental for applica- 
tions, and is thus widely studied in experimental setups. 
Applying an external force on a DES, such as an exter- 
nal magnetic field on a ferromagnetic domain wall or an 
electrical current on a vortex lattice in a type-II super- 
conductor, it is indeed possible to set those systems into 
motion, which is highly non-trivial since they are after all 
emergent structures defined by a shift of the order pa- 
rameter. For example, an external magnetic field favors 
a given magnetization direction, makes the favored ferro- 
magnetic domains grow and thus displace the boundaries 
between the different phases, thus the domain walls effec- 
tively move. 

Ideally we would like to have access to the whole time- 
resolved displacement field u^{t), as in numerical simu- 
lations but also in experiments such as e.g. the imag- 
ing of the imbibition line of a fluid on a disordered sub- 
strate [26], of a crack front along an heterogeneous weak 
plane [27] by ultra-fast CCD camera, or of avalanches 
in ferromagnetic thin films [5]. The quantities of inter- 
est, beyond the displacement field itself u^{t), are usu- 
ally the position of its center of mass 'Ucm(^), its veloc- 
ity v{t) = dtUcuit) and the fluctuations of the velocity 
dtUz{t) itself. However, for a complete characterization, 
the geometrical fluctuations around the center of mass 
Su^{t) = u^{t) — ucuit) should also be addressed, with the 
PDF of the ad hoc dynamical relative displacements at 
fixed time Au^{f) = Su^-^f{t) — Su^{t) and its correspond- 
ing variance B{f^t)^ structure factor S{q^t) or higher mo- 
ments similarly to the statics. 

Analytically, numerically and experimentally there has 
mainly been a focus on the determination of the steady- 
state velocity of the center-of-mass under a constant driv- 
ing force F, with the temperature-dependent velocity-force 
characteristic vt{F) illustrated schematically in Fig. 3. At 
zero temperature, the picture is superficially similar to 
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Figure 3: Velocity- force characteristics: steady-state velocity of the 
center-of-mass position vt as a function of a constant external force 
F. The large arrow shows the displacements of the curves vt{F) 
with an increasing temperature T, coupled with a gradation from 
blue to red curves. 



that of a critical phenomenon, with the steady-state ve- 
locity Vt as an order parameter, the driving force F as the 
controlling parameter and the disorder-dependent critical 
force Fc. Below Fc the DES might relax from a given con- 
figuration with a transient velocity v{t) [28], but then it 
is stuck so vt=o{F < Fc) = 0- Above Fc it acquires a fi- 
nite steady-state velocity vt^o{F > Fc) > through the 
disordered energy landscape, which actually takes off with 
a critical depinning exponent vt=o{F > Fc) ^ {F — Fc)^ . 
With a finite temperature, there are thus three regimes 
of interest for the thermal effects on the dynamics: first 
the depinning regime at F ^ Fc with a possible thermal- 
rounding power law vt{Fc) ^ [29, 30]; secondly the 
large velocity regime at F ^ Fc where vt ^ rjF which 
defines the friction coefficient 77 of the DES, and the pos- 
sible characterization of the nature of this moving phase 
via its dynamical geometrical fluctuations; and finally the 
creep regime at F <C Fc where thermal fluctuations allow 
the system to explore the disordered energy landscape and 
overcome some of its barriers. 

If there was only one characteristic energy barrier A, 
the response of the system to a very small driving force 
would be linear vt{F) ~ e~^^^F ('Thermally Assisted 
Flux Flow'). However such DES actually exhibits a 
disorder-induced metastability and thus there is no fixed 
characteristic energy barrier, and the response is rather a 
stretched exponential: 



vt{F) oc exp 



Uc 

T 



(6) 



This phenomenological 'creep' formula can be related to 
the static quantities assuming that the creepy interface 
moves forward with a succession of avalanches of typi- 
cal size Lc ('Larkin domains'). Assuming that Lc would 
be the characteristic crossover lengthscale in the static 
roughness, and comparing by scaling the elastic, disor- 



der and external force energies at stake at this length- 
scale [31, 32, 33, 15], the creep exponent is predicted to be 
/i = — 2 + 2Ci^M)/(2 — Crm)^ with d the dimension and 
Crm the static asymptotic roughness exponent. This be- 
havior has been spectacularly pointed out for the ID inter- 
face on several order of magnitudes in ferromagnetic DWs 
[4] fi = 0.24 ± 0.04, C = 0.69 ± 0.07, d = 1 and in numeri- 
cal studies [22] . However the content of the typical energy 
Uc and depinning force Fc is still an open question. The 
single-crossover scenario predicts Uc = c^^ jLc the typical 
energy barrier at the lengthscale Lc and Fc = c(^/Ll the 
depinning force of a Larkin domain, with ^ the width of 
the interface, but this prediction has not been reproduced 
yet in numerics even for the ID interface. 

3.3. Methods 

In statics, a given configuration with a displacement 
field Uz is weighted by a Boltzmann weight ex g-'^i^^'^^]/^ 
(with the Boltzmann constant /cb = 1 so that the tempera- 
ture T has actually the units of an energy) , whereas in dy- 
namics the Hamiltonian 1-L [li, V] is replaced by a Martin- 
Siggia-Rose action constructed from a Langevin equation 
for the time-evolution of the displacement field u^{t). The 
statistical averages (O) for the observables O listed in the 
previous sections are then computed analytically with re- 
spect to those statistical weights. 

Pure scaling arguments on the statistical averages of ob- 
servables, via the scaling of the Hamiltonian, of the disor- 
der distribution, of the Langevin equation or of the action, 
can yield very powerful predictions yet they have to be in- 
terpreted carefully, and can often provide a rather a pos- 
teriori short explanation of otherwise painfully computed 
results. For example, the Flory or Imry-Ma 'mean- field' 
argument provides a value Cf quite close to the exact phys- 
ical roughness exponent 

The two main analytical tools used on DES are on one 
hand the Functional Renormalization Group (ERG) where 
the whole disorder correlator (3) evolves under the renor- 
malization procedure [34, 17, 35, 36, 37], and on the other 
hand the Gaussian Variational Method (GVM) as intro- 
duced by Mezard and Parisi on DES, and involving Replica 
to treat the disorder [38, 39]. 

Numerically there are in particular very efficient algo- 
rithm in ID in order to address the dynamics and the 
static of the ID interface, starting from a Langevin equa- 
tion both at zero [40] and at small but finite temperature 
[41, 22]. 

4. Static roughness of a ID interface at finite tem- 
perature 

In this section, we focus on the particular case of the 
ID interface in quenched RB disorder, as an illustration 
of a DES modelling, of the computation and interpreta- 
tion of the static roughness function via GVM and scaling 
arguments, and of the link to the quasistatic creep regime. 
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The study of the ID interface is at the crossroad 
between a fundamental interest in the pecuharities of 
low-dimensional systems, several mappings on related 
statistical-physics problems in the Kardar-Parisi- Zhang 
(KPZ) universality class (including the 1+1 Directed Poly- 
mer (DP)) [42, 43, 44], and experimental realizations of 
effective ID interfaces such as ferroic or ferromagnetic do- 
main walls (DW) in thin films. 

A first example of such experiments are DWs in ferro- 
electric thin films of Pb(Zro.2Tio.8)03 which display out- 
of-plane polarized domains, written and probed by atomic- 
force- microscopy technique [1, 2]. This imaging technique 
probes by definition the surface of the film, which is how- 
ever sufficiently thin (typically 50-60 nanometers) so that 
the measured DWs are effective ID interfaces convoluted 
on the whole thickness of the film. A second example, 
which is believed to be a prototype of the DES model pre- 
sented thereafter, are ferromagnetic DWs in ultrathin films 
of Pt/Co/Pt (just a few atomic layers) with out-of-plane 
magnetization domains [4, 5, 6]. 

The internal structure of those DWs have not been ac- 
cessed, since their typical width (respectively a few unit 
crystal cells and ~ 10 nm) is still below the resolution (re- 
spectively nanometric and micrometric), along with their 
transverse geometrical fluctuations at comparable small 
lengthscales. However, even though their static roughness 
cannot be observed at sufficiently small lengthscales yet, it 
could have indirect consequences on their quasistatic dy- 
namics. Indeed, in order to obtain the phenomenological 
creep formula (cf. section 3.2), the main assumption is 
that there is a single characteristic lengthscale in the stat- 
ics (the Larkin length Lc) which coincides with the typical 
size of avalanches in the quasistatic dynamics. The scal- 
ing prediction for the creep exponent worked surprisingly 
well in ferromagnetic DWs [4] for the roughness and creep 
exponents ( and /i (cf. section 3.2), however the discrepan- 
cies regarding the temperature-dependence of the typical 
energy barrier Uc observed in numerical results [45] sug- 
gested to reexamine the low-temperature static roughness 
at small lengthscales. 

4-1- Full DES model of the ID interface 

Following the DES recipe of the section 2, we consider 
a ID interface {d = m = 1) embedded in a continuous and 
infinite space (z^x) G M^, with implicit ultra-violet and 
infra-red cutoffs ultimately irrelevant in our computations. 
Restricting ourselves to the case without overhangs nor 
bubbles, each configuration is characterized by a univalued 
displacement field with respect to a flat configuration 
defined by the z axis. 

We assume a short-range elasticity in the elastic limit 
as in (1) and a quenched RB weak disorder uncorrelated 
along its internal direction z as in (2). This very particu- 
lar form of elasticity is central for the mappings to other 
statistical physics problems and for the statistical tilt sym- 
metry (STS), a fundamental property for ID interfaces in 
presence of disorder [46, 47, 48]. 




Figure 4: Gaussian profile of a ID interface, described by the nor- 
malized density p^(z, x) as a Gaussian of variance ^ centered at fixed 
z on the interface mean position Uz- 

The disorder correlator (3) is finally taken as a Gaussian 
of variance = 2^^ and amplitude D: 

V47r^ Jr 27r 

So there are two distinct Gaussian hypothesis regarding 
the disorder: on one hand its distribution V \y] is assumed 
to be Gaussian and thus fully described by its disorder cor- 
relator (a physical assumption in the weak disorder limit), 
and on the other hand its disorder correlator itself is cho- 
sen to be a Gaussian. In this choice R{x) decreases suf- 
ficiently fast to be random-bond and it encodes explicitly 
a finite disorder correlation length ^ in a convenient way 
for computations. In an alternative point of view, this 
parameter ^ corresponds to a finite width of the ID inter- 
face, whose spatial extension in the x direction is described 
by a Gaussian density p^{z^x) of variance centered on 
(z^Uz) as illustrated in Fig. 4; if the random-potential is 
5— correlated in both (z, x) directions, the effective disor- 
der Jj^dx • p^{z,x)V{z,x) has precisely the correlator (7). 

4.2. DES characteristic scales and scaling arguments 

In this minimal DES model, there are thus four charac- 
teristic scales: the elastic constant c which gives the energy 
per unit of length of the interface and essentially fixes an 
energy-scale reference; the strength of disorder D or the 
amplitude of the disorder correlator; the finite temperature 
T which quantifies the energy of the surrounding thermal 
bath; and the typical width or disorder correlation length 
^ > or the variance of the disorder correlator. 

The main consequence of the interplay between a finite 
T > and a finite ^ > on the ID interface geometrical 
fluctuations is the existence of two regimes in temperature 
for its static roughness. Those two regimes are separated 
by the single characteristic temperature Tc = {^cDy^^ 
which naturally involves all the scales at stake in the DES 
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model. Physically this can actually be understood by a 
comparison between the microscopic width ^ and an effec- 
tive width ^th due to thermal fluctuations. At high tem- 
perature (T > Tc) the thermal fluctuations erase the exis- 
tence of the microscopic width < ^th), so with the only 
two remaining scales {D^T} one can define a single char- 
acteristic lengthscale r*, presumably the unique crossover 
between two roughness regimes (purely thermal at small- 
lengthscales and dominated by disorder at large length- 
scales). On the contrary, at low temperature (T < Tc) 
^ > 6h so the microscopic width matters and there is a 
whole intermediate roughness regime which leads to the 
definition of at least two characteristic crossover length- 
scales ri and ro, the latter marking the entrance in the 
asymptotic random-manifold regime dominated by disor- 
der. 

The bare scahng of those different quantities r*, ^th, 
Tc and ro can actually be extracted from simple scaling 
arguments, by rescaling the spatial coordinates z = bz and 
X = ax in the whole roughness function: 

5(r;c,AT,0 = J VVV[V]{{u,^r-u,f)^ 

. 2V ■ — ..,.^2.^-H[n,y]/T (8) 

{{Uz^r-U^) )y = 



V - Jp^.g-H[n,y]/T 



Indeed we have to assume first that the random potential 
V 'scales in distribution' consistently with respect to the 
scaling of its disorder correlator, leading thus to 



Since Hei [uz] \ 



-y^el [Uz] 



=1, if we assume that the 
whole Hamiltonian H = Hei + Hdis can be rescaled with 
an overall prefactor, a same scaling for the two parts of 
the Hamiltonian imposes the relation a = D^^'^c~^b^^ with 
the Flory exponent = 3/5 for d = m = 1. Finally we 
can either absorb all the temperature dependence into an 
effective thermal width (fixing the temperature, ^ can still 
be properly neglected so that corresponds to the high-T 
regime): 

B{r; c, D, T, = ^t'h B{r/r,; 1, 1, T = 1, ^/^th) 

T3 7^5 (9) 

WT) = ^, r.(T) = ^ 

or include all the width dependence into a characteris- 
tic temperature (fixing the width, the temperature can be 
pushed arbitrarily low, so that corresponds to the low-T 
regime): 



B{r; c, D, T, = f B{r/ro; 1,1, T/T„ ^ = 1) 



(10) 



Those two raw scalings connect of course, via the condi- 
tions ^th(Tc) = e and r^Tc) = ro(0. 



To give a physical meaning to the different quantities de- 
fined by these scaling arguments, one needs to make addi- 
tional non-trivial assumptions, essentially on the behavior 
of the scaling functions 5(f, 1, 1, 1, f) and 5(f, 1, 1, T, 1), 
e.g. a power-law behavior ~ f'^^. For example the Flory 
exponent could be interpreted as the physical roughness 
exponent ( if there was only one single power-law behavior 
in the roughness, thus the whole system would display a 
true scale invariance at all lengthscales and its Hamilto- 
nian could indeed be rescaled with a scaling factor a ^ b^^ . 
However there are at least two different roughness power- 
law regimes, at small lengthscales a thermal behavior with 
^th(^) = Tr/c ~ (T^/^r-^/^)^ and at large lengthscales the 
asymptotic RM regime with 5asympt(^) ^ (T^r^/^)^. Even 
though the asymptotic roughness exponent of exact value 
Crm = 2/3 can be shown to be robust with respect to 
the addition of a finite width ^ > [49], its temperature- 
dependence (and thorn exponent J)) could a priori be mod- 
ified, as emphasized by the two opposite ERG regimes of 
zero-temperature fixed-point [33] versus high-temperature 
[50]. 

4.3. GVM roughness of the ID interface 

In order to go beyond those scaling predictions and to 
have an explicit roughness function B{r;c,D^T,(_) con- 
necting all the lengthscales at ^ > and T > 0, we have 
computed it in full details in the reference [51] in a Gaus- 
sian Variational Method (GVM), compatible with the 
fairly Gaussian distribution V{Au{r)) observed in numer- 
ics [52]. This scheme has already been applied to DES, 
both in periodic systems [17] as well as manifolds [53], in 
the latter the disorder correlator (7) was assumed to decay 
in a power law, whereas we focused in [51] on the conse- 
quences of its finite variance ~ 

We first used the so-called 'replica trick', well-known 
in the study of spin glasses [39], in order to average first 
over the disorder and so to transform the random part 
Hdis in the full Hamiltonian of one interface (i^i), into an 
effective non-random coupling between n copies of the in- 
terface {u = {ixi, . . . , iin}), in the overall peculiar limit of 
n ^ 0. The effective 'replicated' Hamiltonian H [u] thus 
obtained was then approximated in the statistical thermal 
average via GVM, i.e. it was replaced by a quadratic repli- 
cated Hamiltonian Ho [u] optimized by minimizing its cor- 
responding variational free energy. Note that the GVM ap- 
proximation is the third Gaussian hypothesis in this study 
of the ID interface, based on numerical results [52, 25]. 

In the replica formulation, the Ansatz of the GVM 
self-energy [cr] (li) (u G [O7I]) was taken to be full- 
replica-symmetry-breaking (full-RSB) below a single cut- 
off Vc G [0, 1] and replica-symmetric (RS) above Vc] the 
full-RSB encodes physically the existence of metastability 
at large lengthscales (via a continuum of self-energies) and 
the RS plateau the thermal fluctuations of the interface at 
small lengthscales (via a unique typical self-energy). This 
leads to the following roughness function as a function of 
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10-5 QQQJ T 



can be seen as a directed polymer with one extremity 
fixed at the origin (0,0), growing in 'time' along a tra- 
jectory (t' ^y(t')) G M^, until a time t identified with the 
lengthscale r. With the translations t ^ r and y{t) ^ 
the geometrical fluctuations of the ID interface at a given 
lengthscale r are mapped to the DP's endpoint fluctua- 
tions at a fixed time t, V{Au{r)) ^ V{y{t)) and so their 
variance define the same roughness function: 

B(r) = (Au(r)2) = ^^^^(i) (15) 



Figure 5: Roughness function B(r) at fixed disorder, obtained by 
GVM on the fuh DES model of a ID interface. The large arrow shows 
the displacements of the curves with an increasing temperature T, 
coupled with a gradation from blue to red curves 



In the absence of disorder, the polymer draws a Brown- 
ian random walk and the probability distribution of its 
endpoint is a pure Gaussian given by 



the lengthscale r, illustrated in Fig. 5: 
Bir) 



5k -6 



1 ^ (-f)fe 

^ k=2 

5^7r 1 



Ai(5/6-^c), Ai 



(11) 



+ (l-^c) (12) 



5^7r 
2 X 3^ \T, 



(13) 



(14) 



All the ^-dependence is actually contained in the full- 
RSB cutoff VciO since = {^cD)^/^. At low tempera- 
ture it grows linearly in temperature Vc ^ T/Tc leading to 
the same scaling prediction ro ~ ^^/^c^/^D~'^^^ as in (10), 
whereas at high temperature it saturates at '^c ^ 5/6 and 
we recover ro ^ T^/{cD'^) = r*(T) as in (9). 

Two artefacts of the GVM computation can be directly 
seen in Fig. 5: the asymptotic roughness functions at large 
lengthscales collapse on a single curve 5asympt(^) ~ T^r^/^ 
i.e. Casympt = 3/5 = Cp^ and b = 0. So the GVM pre- 
dicts on one hand the Flory exponent 3/5 instead of the 
exact RM roughness exponent ("rm = 2/3; and on the 
other hand temperature-independent fluctuations at large 
lengthscales in contradiction with the high-T scaling ar- 
gument (9) (except if we assume temperature-dependent 
microscopic parameters {I),c, ^}). 

4-4- Effective DP toy model and its GVM roughness 

An alternative to this GVM computation on a full DES 
model is to use the exact mapping of the ID interface on 
the l-f-l Directed Polymer [54], in order to construct an 
effective 'toy model' containing additional physical ingre- 
dients in a controlled way at a fixed lengthscale^ and only 
then to perform the GVM approximation for the rough- 
ness. 

Assuming as always that there is no overhangs nor bub- 
bles, a segment of length r of the ID interface (z, Uz) G 



exp 



27r 



Tt 



exp 



FthifyY 



(16) 



where the pseudo-free energy Fth(t,y) is the correspond- 
ing term to the elastic Hamiltonian (1). The KB disorder 
explored by the DP translates into an effective free energy 
Fr^{t,y) such that V{y{t))o^ ^-{F,Ut,yHF^{t,y))/T ^ 
fined as an integrated random phase r]{t^ y) in a 'random- 
field' effective formulation: 



Fr^{fy) 



\J y J — oc 



dy-r]{t,y)^cte{t) (17) 



As the corresponding term of the disorder Hamiltonian (2), 
we assumed that at fixed t this effective disorder is Gaus- 
sian, i.e. V [Frj{t)] and V [r]{t)] are fully defined by their 
mean value and their two-point correlator. This is known 
to be exact in the infinite-time limit, at finite temperature 
and zero-width [55]: 



7?(oo, yi)r]{oo, y2) = D • S{yi - ^2) , D = 



cD 



(18) 



[F^(oo,?/i) -F^(oo,?/2)] = D-\yi - y2\ 



where the polymer has completely forgotten its initial con- 
dition and displays a translational-invariance invariance in 
time. So, based on this infinite-time exact result, we as- 
sumed at last that the finite ^ of the microscopic random 
potential V{x^ z) translates simply at finite time t into a 
broadening of the correlator rfrjit)^ of variance 2^^ and 
generic amplitude taken to be Gaussian for the sake of 
computations and equal to (7) with = 2^^. 

Performing again a GVM computation of the roughness, 
with a full-RSB Ansatz for the self-energy below a single 
cutoff Uc G [0, 1] and RS above, we have obtained [51], as 
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Figure 6: Roughness function BDp(t) assuming a fixed effective dis- 
order D = cD/T, obtained by GVM on the Directed Polymer toy 
modeh The large arrows show the displacements of the curves with 
an increasing temperature T, coupled with a gradation from blue to 
red curves. 



illustrated in Fig. 6: 



BMt > tc) 

Bup{t<te) 




^4/3 _ |2 



(19) 



-1/2 



f i e + ^) (20) 



(21) 



predictions to compare with (11)-(14) with the translation 
of 'time' t into the lengthscale r and the identification of 
tc with tq. 

This second GVM procedure is actually carried out at 
fixed 'timeH^ so a priori with 'time'-dependent effective pa- 
rameters and Dt. We can nevertheless assume that first 
the disorder correlation length for the DP is essentially the 
same as for the microscopic random potential i.e. ^ = ^, 
and secondly that D is given by the constant D = cD/T as 
in (18), compatible with the high-temperatures (or zero- 
width) scaling argument (9). This last choice becomes 
however clearly unphysical in the limit T ^ 0, and in fact 
a similar scaling argument on the DP toy model suggests 
that at low temperatures the effective strength of disor- 
der should saturate dX D = cD /Tc in order to recover the 
low- temperatures (or finite- width) predictions (10), both 
for ro(C) and Tc(^). A numerical and analytical study of 
the DP toy model is currently under preparation [49]. 

With the exact 'high-temperature' result D = cD /T ^ 
all the ^-dependence is contained in the full-RSB cut- 
off Uc{C) as in (14), and with the two opposite lim- 
I2 ^ [uc - T/{iD)^/^) and ^ 00 {uc < 3/4) 
recover consistently the high-T scaling predictions 
^Ty{cD^) = n{T) as in (9). 



its 

we 

tc ' 



Figure 7: Log-Log representation of the effective width ^eff of the 
ID interface or DP toy model, as a function of temperature. It is 
defined with respect to the roughness at the beginning of the random- 
manifold asymptotic regime, i.e. at the Larkin length Lc- 



This second GVM computation on the DP toy model ac- 
tually gives the same qualitative predictions as for the full 
ID-interface model (high- versus low-T regimes, interme- 
diate roughness regime at low-T) , but it yields by construc- 
tion of the model the correct asymptotic^roughness expo- 
nent Casympt = 2/3 = and also with D = cD/T di non- 
zero thorn exponent b = -1/3 since 5^'^""^* - 2^-2/3^4/3^ 
compatible with the high-T scaling argument (9). 

4.5. Larkin length and effective width ^eff 

A physical benchmark for the roughness is the Larkin 
length Tc, which is the lengthscale marking the beginning 
of the RM asymptotic regime (Tc ^ r*, ro, tc) and encoded 
in the GVM approach into the maximum value of the self- 
energy [a] {vc) = cr^'^ . This length along the internal di- 
rection of the interface has actually been physically defined 
by Larkin [14] as the lengthscale at which the typical rela- 
tive displacement corresponds to the effective width (,eff of 
the interface: B{Lc) ~ u{Lc)'^ ~ ^^f^. 

Plotting this quantity as a function of temperature in a 
log — log representation (cf. Fig. 7), for both GVM com- 
putations, we can clearly observe the existence of the two 
temperature regimes, by the comparison of the microscopic 
length ^ and a thermal width ^th ^ / {cD) as discussed 
qualitatively with scaling arguments in the section 4.2. 
Those two regimes connect when ^eff(^c) = £,th{Tc) = £, 
defining up to a numerical constant the characteristic 
temperature Tc = (^cT)^/^. This result establishes that 
the two limits ^ ^ and T ^ cannot be innocently ex- 
changed^ as underlined by the existence of the two opposite 
FRG regimes of zero-temperature fixed-point [33] versus 
high-temperature [50]. 

The Larkin length Tc is typically below the resolution of 
the ferromagnetic and ferroelectric DWs experiments men- 
tioned before, and a temperature-dependence study of the 
(asymptotic) roughness and its Jd exponent has not been 
carried out on such ID experimental interfaces yet. Nev- 
ertheless the existence of a low-temperature regime with a 
physically finite width ^ might be relevant even for those 
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experiments, since a crude estimation of Tc for ferromag- 
netic DWs [5] suggests that it could actually be of the order 
of magnitude of room temperature [51]. As for the phe- 
nomenological creep formula discussed previously in the 
section 3.2, the whole scaling argument relating the rough- 
ness and creep exponents and /i, which worked so well for 
ferromagnetic DWs in [4], essentially relies on a compari- 
son of the typical scales at a single characteristic length- 
scale Lc'. a so-called 'Larkin domain' encounters a typical 
energy barrier Uc ^ c^eff/ and it takes a typical minimal 
force Fc ~ c^eff/^c move it. Plugging boldly our scal- 
ing for ^eff(T) and Lc{T) in those expressions, we obtain 
in the high-T regime that Uc ^ T whereas Fc ~ and 
at low-T both those quantities saturate to constants at 
Tc. However, neither of those predictions are compatible 
with the affine behavior for Uc(T) observed numerically 
in [45], but at least now we know that this discrepancy 
might be linked to the appearance at low-T of a whole 
crossover roughness regime at intermediate lengthscales, 
thus jeopardising a priori any argument based on a single 
characteristic lengthscale. 



5. Conclusion 



We have briefly presented in these notes the physics of 
disordered elastic systems. There are many physical real- 
izations of DES, both interfaces or periodic systems. In 
order to probe the disorder-induced glassy properties of 
a DES, one can study simultaneously its statics via its 
geometrical fluctuations (roughness and PDF of the rel- 
ative displacements as a function of the lengthscale) and 
its dynamics via its response to an external force (creep, 
depinning, avalanches, ...). Different characteristic scales 
are at stake, and their competition is encoded at the core 
of the DES formulation - the elastic energy c, the disor- 
der amplitude the disorder correlation length ^ and the 
temperature T - and mirrored in the lengthscales depen- 
dence of observables such as the static roughness or the 
size distribution of avalanches. 

In the low-dimensional case of a ID interface, an addi- 
tional physical ingredient a priori as innocent as the ex- 
istence of a flnite width ^ turns out to be crucial for the 
low-temperature properties in statics and a fortiori also in 
the quasistatic dynamics, as emphasized in the FRG 'zero- 
temperature' flxed point of the renormalization approach. 
We have shown how a combination of scaling arguments 
and concrete GVM computations on a pure static observ- 
able such as the roughness can already yield non-trivial 
elements to a coherent picture of the physics at stake. Ap- 
plication of these ideas in the case of the flnite width in- 
terface to the dynamical quantities is a very challenging, 
and experimentally relevant question, that remains to be 
done. 



Acknowledgements 

This work was supported in part by the Swiss National 
Science Foundation under MaNEP and Division II. 

References 

[1] T. Tybell, P. Paruch, T. Giamarchi, and J.-M. Triscone. Do- 
main wall creep in epitaxial ferroelectric P6(Zro. 2^^0.8)^3 thin 
films. Phys. Rev. Lett. 89, 097601 (2002). 

[2] P. Paruch, T. Giamarchi, and J.-M. Triscone. Domain wall 
roughness in epitaxial ferroelectric PbZro.2Tio.8Os thin films. 
Phys. Rev. Lett. 94, 197601 (2005). 

[3] N. A. Pertsev, D. A. Kiselev, I. K. Bdikin, M. Kosec, and A. L. 
Kholkin. Quasi-one-dimensional domain walls in ferroelectric 
ceramics: Evidence from domain dynamics and wall roughness 
measurements. J. Appl. Phys. 110, 052001 (2011). 

[4] S. Lemerle, J. Ferre, C. Chappert, V. Mathet, T. Giamarchi, 
and P. Le Doussal. Domain wall creep in an Ising ultrathin 
magnetic film. Phys. Rev. Lett. 80, 849 (1998). 

[5] V. Repain, M. Bauer, J. -P. Jamet, J. Ferre, A. Mougin, C. 
Chappert, and H. Bernas. Creep motion of a magnetic wall: 
Avalanche size divergence. Europhys. Lett. 68, 460 (2004). 

[6] P. J. Metaxas, J. P. Jamet, A. Mougin, M. Cormier, J. Ferre, 
V. Baltz, B. Rodmacq, B. Dieny, and R. L. Stamps. Creep 
and fiow regimes of magnetic domain-wall motion in ultrathin 
Pt/Co/Pt films with perpendicular anisotropy. Phys. Rev. Lett. 
99, 217208 (2007). 

[7] S. Moulinet, C. Guthmann, and E. Rolley. Roughness and dy- 
namics of a contact line of a viscous fiuid on a disordered sub- 
strate. Eur. Phys. J. E 8, 437 (2002). 

[8] S. Santucci, K. J. Mal0y, A. Delaplace, J. Mathiesen, A. Hansen, 
J. O. Haavig Bakke, J. Schmittbuhl, L. Vanel, and P. Ray. 
Statistics of fracture surfaces. Phys. Rev. E 76, 016104 (2007). 

[9] G. Blatter, M. V. Feigel'man, V. B. Geshkenbein, A. I. Larkin, 
and V. M. Vinokur. Vortices in high-temperature superconduc- 
tors. Rev. Mod. Phys. 66, 1125 (1994). 
[10] G. Coupler, C. Guthmann, Y. Noat, and M. S. Jean. Local 
symmetries and order-disorder transitions in small macroscopic 
wigner islands. Phys. Rev. E 71, 046105 (2005). 
[11] T. Giamarchi. Electronic glasses. arXiv:cond-mat/0403531 
(2004). 

[12] G. Griiner. The dynamics of charge-density waves. Rev. Mod. 

Phys. 60, 1129 (1988). 
[13] S. Brazovskii and T. Nattermann. Pinning and sliding of driven 

elastic systems: from domain walls to charge density waves. 

Adv. Phys. 53, 177 (2004). 
[14] A. I. Larkin. Effect of inhomogeneities on the structure of 

the mixed state of superconductors. Sov. Phys. JETP 31, 784 

(1970). 

[15] T. Giamarchi, A. B. Kolton, and A. Rosso. Dynamics of disor- 
dered elastic systems. Lecture Notes in Physics 688, 91 (2006); 
arXiv:cond-mat/0503437 (2005). 

[16] T. Giamarchi. Encyclopedia of Complexity and Systems Sci- 
ence, chapter Disordered Elastic Media, p. 2019, Springer 
(2009). 

[17] T. Giamarchi and P. Le Doussal. Elastic theory of fiux lattices 
in the presence of weak disorder. Phys. Rev. B 52, 1242 (1995). 

[18] K.-J. Kim, J.-C. Lee, S.-M. Ahn, K.-S. Lee, C.-W. Lee, Y. J. 
Cho, S. Seo, K.-H. Shin, S.-B. Choe, and H.-W. Lee. Interdi- 
mensional universality of dynamic interfaces. Nature 458 , 740 
(2009). 

[19] J. F. Joanny and P. G. de Gennes. A model for contact angle 

hysteresis. J. Chem. Phys. 81, 552 (1984). 
[20] P. L. Doussal and K. J. Wiese. Elasticity of a contact-line 

and avalanche-size distribution at depinning. Phys. Rev. E 82, 

011108 (2009). 

[21] A. Rosso and W. Krauth. Origin of the roughness exponent in 
elastic strings at the depinning threshold. Phys. Rev. Lett. 87, 
187002 (2001). 



10 



[22] A. B. Kolton, A. Rosso, T. Giamarchi, and W. Krauth. Creep 

dynamics of elastic manifolds via exact transition pathways. 

Phys. Rev. B 79, 184207 (2009). 
[23] W. T. M. Irvine, V. Vitelli, and P. M. Chaikin. Pleats in crystals 

on curved surfaces. Nature 468, 947 (2010). 
[24] V. Lecomte, S. E. Barnes, J. -P. Eckmann, and T. Giamarchi. 

Depinning of domain walls with an internal degree of freedom. 

Phys. Rev. B 80, 054413 (2009). 
[25] A. Rosso, R. Santachiara, and W. Krauth. Geometry of gaussian 

signals. J. Stat. Mech. 2005, L08001 (2005). 
[26] S. Santucci, R. Planet, K. J. Mal0y, and J. Ortm. Avalanches 

of imbibition fronts: Towards critical pinning. Europhys. Lett. 

94, 46005 (2011). 
[27] K. T. Tallakstad, R. Toussaint, S. Santucci, J. Schmittbuhl, and 

K. J. Mal0y. Local dynamics of a randomly pinned crack front 

during creep and forced propagation: An experimental study. 

Phys. Rev. ESS, 046108 (2011). 
[28] A. B. Kolton, A. Rosso, E. V. Albano, and T. Giamarchi. Short- 
time relaxation of a driven elastic string in a random medium. 

Phys. Rev. B 74, 140201 (2006). 
[29] S. Bustingorry, A. B. Kolton, and T. Giamarchi. Thermal 

rounding of the depinning transition. Europhys. Lett. 81, 26005 

(2008). 

[30] S. Bustingorry, A. Kolton, A. Rosso, W. Krauth, and T. Gi- 
amarchi. Thermal effects in the dynamics of disordered elas- 
tic systems. Physica B: Condensed Matter 404, 444 (2009); 
arXiv:0809.4671vl. Proceedings of the International Workshop 
on Electronic Crystals on Strongly Correlated Electron Systems. 

[31] T. Nattermann. Interface roughening in systems with quenched 
random impurities. Europhys. Lett. 4, 1241 (1987). 

[32] L. B. loffe and V. M. Vinokur. Dynamics of interfaces and 
dislocations in disordered media. J. Phys. C 20, 6149 (1987). 

[33] P. Chauve, T. Giamarchi, and P. Le Doussal. Creep and depin- 
ning in disordered media. Phys. Rev. B 62, 6241 (2000). 

[34] D. S. Fisher. Interface fluctuations in disordered systems: 5 — e 
expansion and failure of dimensional reduction. Phys. Rev. Lett. 
56, 1964 (1986). 

[35] P. Le Doussal, K. J. Wiese, and P. Chauve. Functional renor- 
malization group and the fleld theory of disordered elastic sys- 
tems. Phys. Rev. EQ9, 026112 (2004). 

[36] K. J. Wiese and P. L. Doussal. Functional renormalization for 
disordered systems, basic recipes and gourmet dishes. Markov 
Processes Relat. Fields 13, 777 (2007). 

[37] P. L. Doussal. Exact results and open questions in flrst principle 
functional RG. Ann. Phys. 325, 49 (2010). 

[38] T. Castellani and A. Cavagna. Spin-glass theory for pedestrians. 
Journal of Statistical Mechanics: Theory and Experiment page 
P05012 (2005). 

[39] M. Mezard, G. Parisi, and M. A. Virasoro. Spin Glass The- 
ory and Beyond volume 9 of World Scientific Lecture Notes in 
Physics. World Scientiflc (1987). 

[40] A. Rosso and W. Krauth. Roughness at the depinning threshold 
for a long-range elastic string. Phys. Rev. EQh, 025101 (2002). 

[41] A. B. Kolton, A. Rosso, T. Giamarchi, and W. Krauth. Dynam- 
ics below the depinning threshold in disordered elastic systems. 
Phys. Rev. Lett. 97, 057001 (2006). 

[42] 1. Corwin. The Kardar-Parisi-Zhang equation and universality 
class. arXiv:1106.1596 (2011). 

[43] J. Krug. Origins of scale invariance in growth processes. Adv. 
Phys. 46, 139 (1997). 

[44] T. Halpin-Healy and Y.-C. Zhang. Kinetic roughening phenom- 
ena, stochastic growth, directed polymers and all that. Aspects 
of multidisciplinary statistical mechanics. Phys. Rep. 254, 215 
(1995). 

[45] A. B. Kolton, A. Rosso, and T. Giamarchi. Creep motion of 
an elastic string in a random potential. Phys. Rev. Lett. 94, 
047002 (2005). 

[46] E. B. U. Schulz, J. Villain and H. Orland. Thermal fluctuations 
in some random fleld models. J. Stat. Phys. 51, 1 (1988). 

[47] T. Hwa and D. S. Fisher. Anomalous fluctuations of directed 
polymers in random media. Phys. Rev. B 49, 3136 (1994). 



[48] D. S. Fisher and D. A. Huse. Directed paths in a random po- 
tential. Phys. Rev. B 43, 10728 (1991). 

[49] E. Agoritsas, V. Lecomte, and T. Giamarchi (in preparation). 

[50] S. Bustingorry, P. Le Doussal, and A. Rosso. Universal high- 
temperature regime of pinned elastic objects. Phys. Rev. B 82, 
140201 (2010). 

[51] E. Agoritsas, V. Lecomte, and T. Giamarchi. Temperature- 
induced crossovers in the static roughness of a one-dimensional 
interface. Phys. Rev. B 82, 184207 (2010). 

[52] T. Halpin-Healy. Directed polymers in random media: Proba- 
bility distributions. Phys. Rev. A 44, R3415 (1991). 

[53] M. Mezard and G. Parisi. Replica fleld theory for random man- 
ifolds. J. Phys. II, 809 (1991). 

[54] M. Mezard and G. Parisi. Manifolds in random media: two 
extreme cases. J. Phys. 12, 2231 (1992). 

[55] D. A. Huse, C. L. Henley, and D. S. Fisher. Huse, Henley, and 
Fisher respond. Phys. Rev. Lett. 55, 2924 (1985). 



Contents 



1 Introduction 1 

2 DES modelling: a recipe 1 

2.1 Dimensionality and class of DES 1 

2.2 Physical space of coordinates (z, x) 2 

2.3 Univalued displacement field 2 

2.4 Elasticity 2 

2.5 Disorder 2 

2.6 Internal structure of DES 3 

3 Observables as probe of disorder: statics ver- 
sus dynamics 3 

3.1 Statics: geometrical fluctuations and 
roughness 3 

3.2 Dynamics: velocity-force characteristic ... 4 

3.3 Methods 5 

4 Static roughness of a ID interface at finite 
temperature 5 

4.1 Fuh DES model of the ID interface 6 

4.2 DES characteristic scales and scaling argu- 
ments 6 

4.3 GVM roughness of the ID interface 7 

4.4 Effective DP toy model and its GVM rough- 
ness 8 

4.5 Larkin length and effective width ^gff • • • • 9 

5 Conclusion 10 



11 



